The Topological Relation Between Bulk Gap Nodes and Surface Bound States 
Application to Iron-based Superconductors 
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In the past few years materials with protected gapless surface (edge) states have risen to the central 
stage of condensed matter physics. Almost all discussions centered around topological insulators and 
superconductors, which possess full quasiparticle gaps in the bulk. In this paper we argue systems 
with topological stable bulk nodes offer another class of materials with robust gapless surface states. 
Moreover the location of the bulk nodes determines the Miller index of the surfaces that show (or not 
show) such states. Measuring the spectroscopic signature of these zero modes allows a phase-sensitive 
determination of the nodal structures of unconventional superconductors when other phase-sensitive 
techniques are not applicable. We apply this idea to gapless iron-based superconductors and show 
how to distinguish accidental from symmetry dictated nodes. We shall argue the same idea leads 
to a method for detecting a class of the elusive spin liquids. 
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The interest in topological insulators and superconduc- 
tors lie mainly in studying their robust gapless surface 
states 1,2 . Systems with bulk gap nodes are usually ne- 
glected for it is generally felt bulk nodes render surface 
states ill defined. Gap nodes in the form of band cross- 
ing are singularities in the electronic structure. Examples 
include the Dirac nodes in graphene, the d-wave super- 
conducting gap nodes of cuprates and the Weyl nodes 3,4 . 
These singularities have dramatic effects on the bound- 
ary electronic structure. For example, the zigzag edge of 
graphene has zero energy (E=0) flat bands 5,6 , and the 
{110} surfaces of the cuprate superconductors have zero 
bias Andreev bound states (ZBABS) 7,8 . Recently it has 
been shown that the Weyl nodes will give rise to surface 
"Fermi arcs" 4 . We shall argue the above examples are 
special cases of a more general fact, namely translational 
invariant systems with topologically stable bulk nodes also 
possess protected gapless surface bound states. In the 
presence of disorder the otherwise degenerate zero en- 
ergy band broadens into finite width, but does not lose 
the signature of gapless surface bound state. 

An important class of "nodal material" is nodal super- 
conductor. Usually, the presence of gap nodes reflects the 
overall repulsive nature of the (effective) Cooper pairing 
Hamiltonian. The location of the gap nodes provide de- 
tailed information about the momentum dependence of 
the pairing interaction. This information is necessary for 
the understanding of the pairing mechanism of nodal, 
hence unconventional, superconductors. Among these 
superconductors the pairing symmetry of iron based su- 
perconductors (FeSCs) is still very much under debate. 
Experimentally roughly half of the materials in this fam- 
ily show the evidence of gap nodes. Examples include 
LaFePO 9 , BaFe 2 (Asi_ x P :E ) 2 10 - :L3 , KFe 2 As 2 14,15 , LiFeP 16 
etc. The nodal structures of the majority in this group 
are unknown. The notable exception is BaFe 2 (Asi_ x P x ) 2 
whose node locations have been determined by angle- 
resolved photoemission spectroscopy 13 . Because of the 
complex fermiology and the likelihood of "accidental 



nodes" these materials are not accessible to conventional 
phase-sensitive probes. Finding experimental means to 
pin down the location of the nodes for these materials has 
been a major experimental challenge. In the following we 
propose to study the correlation between the orientation 
of the surfaces with the presence/absence of zero energy 
peaks in surface sensitive spectroscopy. We show this 
method does not suffer the above difficulty. We will il- 
lustrate the idea by focusing on LaFePO. Near the end 
we shall briefly discuss gapless non-centrosymmetric su- 
perconductors, and generalize the idea presented here for 
the detection of a class of spin liquids. 

Bulk nodes and zero energy surface bound 
states. The gap nodes of a d-dimensional quasi- 
particle Hamiltonian can have non-zero dimension 
q. For examples q = for point node and q = 1 
for line node etc. A prerequisite of a stable nodal 
structure is in d — q — 1 dimension the symmetry of the 
Hamiltonian 17,18 at a generic k point admits topological 
insulator/superconductors 19,20 . If a nodal structure is 
stable, and if d — q — 1 > 0, then suitably oriented 
boundaries of this material will have E=0 bound state. 
Generically the dimension of the sub boundary BZ pos- 
sessing E=0 bound states is dE=o = q + 1. For graphene 
d = 2, q = 0; for cuprates d = 3, q = 1; and for the Weyl 
semimetal 4 d = 3, q = 0. For the purpose of this paper 
we focus on the cases where the band topology of the 
d — q — 1 dimensional system is classified by the group of 
integers 17 . Under this condition a winding number 17 (or 
vorticity) can be assigned to each independent nodal 
manifold (e.g. point or line). As long as the projection 
of opposite-winding-number nodal manifolds does not 
completely overlap in the boundary Brillouin zone (BZ), 
there will be gapless surface bound states. Special cases 
of the above argument have already been put forward 
in Ref. 4 ' 21,22 , especially in Volovik et al. 23 . A more 
illustrative example is given in Appendix A. 

Application to LaFePO. LaFePO is a stoichiomet- 
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FIG. 1: Caricatures of the Fermi surfaces of LaFePO and four 
examples of nodal structures. The blue/red curves represent 
the hole/electron Fermi surfaces. Red dots with ± signs indi- 
cate nodes with ± sign of winding numbers. (a,b) Two differ- 
ent s-wave gap nodes (red dots), (c) d x 2_ y 2 nodes and (d) p x 
nodes. The black circle in part (a) is a close loop enclosing a 
gap nodes. 



ric compound, it is clean enough for quantum oscilla- 
tion to be observed 30 . According to Ref. 30 LaFePO has 
nearly two dimensional Fermi surfaces: two hole pock- 
ets at the center and two electron pockets on the face 
of the unfolded BZ (for a caricature see Fig. la). In 
the following we shall used 2D notations. LaFePO is 
the first FeSC where strong evidence of line nodes have 
been reported 9 . At the present time there is no direct 
information about the location of the line nodes. In ad- 
dition, early Knight shift measurement has revealed a 
rather puzzling result. Unlike conventional singlet-paired 
superconductors, the Knight shift increases rather than 
decreases below T c 31 ' 32 . This unusual behavior calls for 
the examination of possible nodal triplet pairing. 

Theoretically it has been proposed that the singlet 
gap function of LaFePO possesses the so-called "acciden- 
tal nodes" 33-35 where the nodal gap function transforms 
without any sign change under all point group operations. 
A schematic illustration of the proposed nodal structure 
is shown in Fig. la. First let's discuss the stability of 
these accidental nodes. Consider a loop encircling one of 
nodes (e.g. the black circle in Fig. la). Like the cuprates 
the BdG Hamiltonian defined on the loop is character- 
ized by an even integer winding number, and has a fully 
gapped spectrum. Thus the loop has a non-trivial band 
topology. If the node is removed from the interior of the 
circle the winding number must change, which requires 
the spectrum on the circle to become gapless, i.e. the 



node must pass through the circle. Of course the above 
statement applies to every loop enclosing the node. This 
leads to the "vorticity conservation law" - a single node 
can be displaced but not annihilated. Annihilation can 
occur only when two opposite vorticity nodes hit each 
other. Therefore although "accidental" these nodes are 
stable against small perturbations. 

Because the opposite vorticity nodes distribute sym- 
metrically around the (10), (01) and (11), (11) axes in 
Fig. la, there will be no ZBABS on the edges. This is in 
sharp contrast to the d x 2_ y 2 nodes (Fig. lc) which pos- 
sess ZBABS on the (11) and (11) edges, as well as the 
d xy nodes which have ZBABS on the (10) and (01) edges. 
The lowest index edge that does show ZBABS are the 
{21} edges. (The corresponding orientation of the real 
surfaces can be obtained by adding an extra zero, e.g., 
(11) ->■ (110) etc.) Using the tight-binding fit of LDA 
bandstructure adapted from Ref. 33 ' 34 and the following 
BdG Hamiltonian 

#LaFcPo(k) = e(k) <g> cr (g> t 3 + A(k)/ 5x5 ® a <g> n, (1) 

we calculated the surface local density of states for a 
number of different singlet pairing symmetries and sur- 
face orientations. In Eq. (1) e(k) is the 5x5 bandstruc- 
ture matrix, A(k) = 2Ao(cosfca; +cos/c y ) and /s X 5 is the 
identity matrix in the orbital space. (The calculation 
uses the iterative Green's function method of Ref. 39 .) A 
few examples of the results for nodal singlet gap function 
are shown in Fig. 2a-d. 

Of course the s-wave accidental nodes can fall on the 
central hole pockets too, an example is shown in Fig. lb. 
(Not dictated by symmetry, there is no reason for ac- 
cidental nodes to occur on both hole Fermi surfaces.) 
Due to the reflection symmetry about the {10} and {11} 
axes, the projection of the opposite vorticity nodes coin- 
cide, hence there will be no ZBABS on those surfaces. In 
contrast there is no reflection symmetry about the {21} 
directions hence generically these edges will have ZBABS 
. This argument holds when the accidental nodes exist on 
both hole pockets and/or on the electron pockets. 

Because the unusual temperature dependence of the 
Knight shift we have also considered triplet pairing. The 
BdG Hamiltonian for this case is given by 

#LaFcP ( k ) = e( k ) ® 0-0 <8 T 3 + / 5x5 ® A(k) ® Tl, (2) 

where A(k) = d(k) ■ a with d(k) = — d(— k). In general 
if there is no further restriction on d(k) the node should 
not be topologically stable. This is because for a loop 
consists of generic momenta around a node, the mapping 
from k to (e(k),d(k)) (e(k) 2 + \d{k)\ 2 ^ 0) has the homo- 
topy group iri(S 3 ) = 0. The situation can change if we 
impose further constraints on the form of d(k). For ex- 
ample, ignoring the spin-orbit interaction, if we require 
d(k) = /(k)rfo where /(k) is a partner function of the 
point group irreducible representation, the node becomes 
stable, and is characterized by the winding number of the 
AIII class 17 . Under such condition we expect boundary 
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FIG. 2: The surface local density of states N(lj) for some 
of the pairing symmetries in Table I. For each symmetry 
two surface orientations are shown: one with zero bias peak 
and the other without. The unit of the density of states is 
arbitrary. The Ao used in this calculation is O.OleV. Thermal 
broadening with temperature T = 0.4meV is used. 



zero energy states. For simplicity we considered only the 
following p-wave pairings with <i(k) = 2Ao sin k x z (p x 
pairing), 2Arjsinfc a z (p y pairing) and 2Ao sin(/c x ±k y )z 
(Px± y pairing). 

In Table I we summarize the results for eight different 
time reversal symmetric nodal pairing scenarios. In addi- 
tion to the four triplet gap functions discussed above we 
have also studied nodal s, A(k) = 2A (cosfc :E + cosfc !y )ero; 
d x 2_ y 2, A(k) = 2Ao(cosfca; — cosk y )ao and d xy , A(k) = 
4Aq sin k x sin k y <Jo (here Ao = O.OleV). Moreover we 
have also considered the horizontal ring nodes observed 
in Ref. 13 for BaFe2(Asi_ 2; P a ;)2- Surface orientation wise 
we considered (010), (100), (110), (110) and (210). Since 
there are no identical rows, by studying the dependence 
of the zero-bias peak on the surface orientation one can 
differentiate different types of nodal pairing symmetries. 

Although our discussion is in the context of LaFePO, 
our results should apply to other iron-based superconduc- 
tors with the same lattice symmetry and time-reversal 
symmetry as well. For example the proposed d x i_ y i- 
wave pairing in KFe2As2 36 should have the same surface 
orientation dependence of ZBABS as the cuprates. 

Effects of disorder. A natural question one might 
ask in reading this work is how protected these E=0 sur- 
face states really are. Because r 2 , the second Pauli ma- 
trix in the Nambu space, anticommutes with all Hamil- 
tonians we have discussed, the E=0 states can be made 
eigenstates of it. It can be shown that for a fixed winding 
number the E=0 surface states with different ti eigenval- 



TABLE I: V" denotes the presence of ZBABS and "x" the 
absence of it. The pairing symmetries are defined in main 
text. 
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ues localize at the opposite ends of the sample. Revers- 
ing the winding number exchanges these localized states. 
Thus local disorder potential can only mix E=0 states 
at surface momenta where opposite winding numbers oc- 
cur. It turns out time reversal symmetry requires the ID 
winding number to change sign from surface momentum 
k to — k. As the result there are equal, opposite, inter- 
vals in the boundary BZ with opposite winding numbers. 
Smooth surface potentials do not mix them but sharp 
impurity potentials do. In the presence of the mixing 
the delta function density of states broadens into a peak 
centered around E = 0. The width of the peak is pro- 
portional to the disorder strength. Therefore although 
the singular density of states at E=0 is not protected, 
a broadened version of it is. Interestingly the p x ^ y -w&ve 
pairing studied in conjunction with LaFePO evades dis- 
order broadening hence is strictly protected. This is be- 
cause opposite winding number occurs for opposite spin. 
As long as the disorder potential does not flip spins the 
E = singularity is unaffected. Of course the ability 
to see the zero bias conductance peak 24 ' 25 and its sur- 
face orientation dependence 26-29 in real experiment on 
cuprates, where disorder should mix ZBABS, attests for 
the robustness of the phenomenon discussed here. 

Nodal non-centrosymmetric superconductors. 
Another interesting application of the idea presented 
here is to nodal non-centrosymmetric superconductors 37 . 
In these systems due to the spin-orbit interaction the 
Fermi surfaces are spin non-degenerate. This is shown 
by the black and gray curves in Fig. 3a. Several known 
non-centrosymmetric superconductors exist near an an- 
tiferromagnetic phase. These include CePtsSi, CeRhSi3, 
CeIrSi3 and CeCoGe3 37 . For these systems it is natu- 
ral to expect the superconducting gap function to pos- 
sess nodes (in particular d-wave nodes). This is indeed 
substantiated by experimental findings 38 . Therefore ac- 
cording to earlier arguments there should be surface E=0 
bound states. We consider the following simple model 
BdG Hamiltonian 

H(k) =e(k)ao ® r 3 - A<y(k) ■ a ® t 3 

+ Aq(cos k x — cos k y )[aao + /?7(k) - a] ® T\. 



Here e(k) 



•cosfea;— cosky — /i, 7*(k) = (sin k y . — sin k x ) 
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FIG. 3: (a) The spin non-degenerate Fermi surfaces (the black 
and gray curves) of the non-centrosymmetric superconductor 
defined in Eq. (3). The nodes are designated by the red dots, 
(b) The bandstructure of the {11} edges. The edge states 
exist between the momenta marked by the yellow arrows. The 
parameters we used in producing Fig. 3b are A = 0.3, fi = 
0.45, A = 0.1 and a = 1, P = 0. 



and a = ( (71,(72). The bandstructure on the {11} edges 
of this model Hamiltonian is shown in Fig. 3b. There 
are E=0 flat bands in different portions of the edge BZ. 
Generally, when exist, such flat bands are two fold de- 
generate, i.e., on each edge there are two E = bound 
states. The only exception are the flat bands (marked 
by the red thick line segments in Fig. 3b) in the momen- 
tum intervals between the yellow arrows. The momenta 
of these two intervals are opposite to each other (up to 
reciprocal lattice vector). In addition the zero energy 
bound states there are non- degenerate. The E=0 quasi- 
particle operator associated with the two different red in- 
tervals are hermitian conjugate of each other. In Rcf. 22 
similar surface flat bands are found in the nodal states 
bordering non-centrosymmetric (fully gapped) topologi- 
cal superconductors. 

Detecting nodal spin liquids. The search for spin 
liquid is a focused interest in current condensed matter 
physics. The difficulty with finding a spin liquid is pre- 
cisely why they are novel - they are "featureless" in con- 
ventional sense. Among different spin liquids there is a 
class that are nodal (to ensure their stability we focus on 



the ones having Z2 gauge symmetry). For example such 
a spin liquid was proposed for the pseudogap state of 
the cuprates 40,41 . Like nodal superconductors described 
above, nodal spin liquid should have surface "spinon" 
flat bands. If exist these flat bands will give rise to Curie 
behavior of local magnetic susceptibility hence allows for 
experimental detection, for example the low temperature 
susceptibility will diverge as 1/T. Of course, as discussed 
earlier, disorder broadens the singular density of state at 
E = into a peak hence renders the local magnetic sus- 
ceptibility finite. In addition one might need to worry 
about the effect of residual spinon interaction on these 
flat bands. We believe if the residual spinon interaction 
is repulsive surface ferromagnetism will result. 

Discussions. Experimentally one can prepare dif- 
ferently oriented surfaces by thin film growth and per- 
form scanning tunneling microscopy (STM) measure- 
ments to measure the local density of states associated 
with ZBABS. The advantage of using STM is it can probe 
~ nm length scale. Therefore even when the surface is 
not ideal, it can pick up signals from surfaces with dif- 
ferent local orientations. Thus by scanning the tip it is 
possible to acquire the required data from a single sam- 
ple. 

After completion of this work we noticed a very re- 
cent experiment observing zero-bias conductance peak 
for FeSe . 3 Teo.7 42 - 
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Appendix A: The criterion illustrated by the 
cuprate example 

The idea behind the criterion presented in the main 
text is best illustrated by using the cuprate supercon- 
ductor as an example. The Bogoliubov-de Gennes (BdG) 
Hamiltonian of the cuprate superconductor read 

#cupr a te S (k) = e(k)cr <g> T 3 + A(k)cr ® 71, (Al) 

where Co is the identity matrix acting in the spin space 
and Ti ; 3 are 2x2 Pauli matrices in the Nambu space, e(k) 
is the normal state dispersion satisfying e(— k) = e(k) and 
A(k) is the d-wave gap function. (Since the cuprates 
are quasi two dimensional materials, we shall use two 
dimensional notations in the following discussions.) The 
Fermi surface and the gap nodes are shown in Fig. 4a, 
therefore d = 2, q = 0. In the same figure, the normalized 
vector 

n(k) = (e(k),A(k))/v/e(k) 2 + A(k) 2 (A2) 
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FIG. 4: (a) The 2D nodal band structure and its projections 
on (11) and (01) surfaces for a model d-wave superconductor. 
The first BZ and one extended zone are drawn. The solid 
black curves denote the Fermi surface and the shaded region 
is filled in the normal state. Red arrows are the unit vectors 
in Eq. (A2). The blue dots with ± signs represent nodes with 
vorticity ±2 respectively. Top left: the slanted thin black 
line segment is the surface BZ of the (11) edge. The black 
arrow with letter "P" indicates the direction of projection. 
Thick red line segments mark the surface momenta with zero 
energy bound states. Top right: the thin black horizontal line 
segment is the BZ of the (01) edge. The ± vorticity nodes 
overlap after projection. The vorticity of the node enclosed by 
the parallelogram is equal to the winding number difference 
along the two side blue segments because the winding along 
the top and bottom segments cancel due to the periodicity 
in momentum space. This can be explicitly seen by following 
the turning of the red arrows (the actual winding number is 
twice of the winding shown by the arrows, due to the spin 
degeneracy), (b) The (11) edge bandstructure. The surface 
flat bands are marked red. In constructing this figure we have 
used e(k) = — cosk x — cosk y — /j, and A(k) = Ao(cosfca; — 
cosfcy) in Eq. (Al). Here fi = 0.45, A = 0.1. 



is plotted as a function of k over two BZs (see red ar- 
rows). Inspecting these arrows one notices each node is a 
"vortex" in n(k). Around each vortex the arrows exhibit 
non-zero winding. The total winding number associated 
with each node is given by 



w = — <b dp • [m(k)V k n 2 (k) - n 2 (k)V k m(k)]. (A3) 



(The extra factor of 2 is due to spin degeneracy) . Clearly 
each node is characterized by an even integer wind- 
ing number. The BdG Hamiltonian defined on all one 
(= d — q — 1) dimensional loop enclosing the node are 
topologically nontrivial. 



Now consider the bandstructure projected along the 
(11) direction. For each transverse momentum k along 
(11) we have a ID chain running in (11). So long as k 
does not coincide with the projection of the nodes the 
spectrum is fully gapped and characterized by the wind- 
ing number defined in Eq. (A3). Any two chains whose 
k straddle the projection of a node their winding num- 
bers must differ by ±2 (see Fig. 4a captions), hence at 
least one of them is topologically non-trivial and possess 
E = end states when the boundary condition along (ll) 
changes from closed to open. This implies E=0 bound 
states exists for intervals of k. Therefore dE=o is indeed 
q + 1 = 1. An example of the (11) boundary bandstruc- 
ture is shown in Fig. 4b. The k intervals showing the 
flat bands are represented by the thick red line segments 
in the top left corner of Fig. 4a. The only edges which 
do not possess ZBABS are the {10} (Miller's notation is 
used) edges where the projection of positive and negative 
nodes overlap (see top right corner of Fig. 4a). For the 
real material d = 3, q = 1 and the only modification is 
dE=o changes from 1 to 2. 
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